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Introduction 


Initial  studies  of  the  effects  of  plasma  density  gradients  on  the  diffu¬ 
sion  dominated  regions  of  the  ionosphere  concentrated  on  the  application  of 
plasma  transport  equations  to  solve  for  plasma  density  altitude  profiles  for 
the  case  of  the  diffusive  equilibrium  distribution  of  multi-ion  species  (cf 
Mange,  1960).  In  this  case,  the  pressure  gradient  existing  between  the 
F-region  peak  and  high  altitudes  results  in  the  ions  distributing  themselves 
such  that  the  polarization  electric  field  set  up  by  the  pressure  gradient  is 
just  balanced  by  the  gravitational  force  on  the  ions.  Under  these  circum¬ 
stances  the  outward  flux  of  ions  is  zero.  When  the  concept  of  open  or  extend¬ 
ed  geomagnetic  field  lines  became  accepted,  it  was  realized  that  the  possibil¬ 
ity  existed  for  ions  to  continuously  flow  from  the  earth  into  the  distant 
regions  of  the  magnetosphere  linked  by  geomagnetic  field  lines  to  the  polar 
ionosphere.  This  condition  is  one  of  the  solutions  of  the  ion  density  distri¬ 
bution  under  diffusion  dominated  conditions  and  is  known  as  "dynamic  equili¬ 
brium"  as  opposed  to  "diffusive  equilibrium"  discussed  earlier.  Both  types  of 
equilibrium  represent  steady  state  solutions  to  the  plasma  transport  equations 
resulting  in  the  requirement  that  the  outward  flux  is  a  constant,  zero  in  the 
case  of  diffusive  equilibrium  and  some  finite  value  in  the  case  of  dynamic 
equilibrium  (Banks  and  Kockarts  1973;  Bauer,  1973).  Boundary  conditions  on 
the  rate  of  production  of  the  outflowing  ions  and  Coulomb  collisions  with  sta¬ 
tionary  ions  result  in  an  upper  limit  to  the  outward  flux.  Numerous  theoreti¬ 
cal  models  have  been  developed  over  the  last  decade  to  describe  the  dynamic 
equilibrium  condition;  these  models  include  hydrodynamic  models  (Banks  and 
Holzer,  1968,  1969a,  1969b;  Marubashi ,  1970;  Banks,  1973;  Bailey  and  Moffett, 
1974;  Banks  et  al,  1974,  1976;  Strobel  and  Walker,  1972;  Raitt  et  al,  1973, 


1977,  1978a,  1978b;  Schunk  et  al,  1978;  Ottley  and  Schunk,  1980),  hydromagnet- 
ic  models  (Holzer  et  al,  1971),  kinetic  models  (Lemaire  and  Scherer,  1970, 
1971,  1972,  1973;  Lemaire,  1972)  and  models  based  on  generalized  transport 
equations  (Schunk  and  Watkins,  1979,  1981,  1982). 

This  strong  emphasis  on  the  description  of  the  dynamic  equilibrium  cases 
has  resulted  in  the  diffusive  equilibrium  case  being  somewhat  neglected  since 
the  earlier  work  referenced  above.  In  particular  the  extension  of  model  cal¬ 
culations  of  both  the  dynamic  equilibrium  case  and  the  diffusive  equilibrium 
case  along  a  magnetic  flux  tube  all  the  way  from  the  ionosphere  to  the  equato¬ 
rial  plane  has  been  neglected.  Chiu  et  al  (1979)  corrected  this  deficiency 
for  the  diffusive  equilibrium  case  by  describing  a  model  for  plasma  density 
distribution  from  the  ionosphere  to  the  equatorial  plane  allowing  for  magnetic 
flux  tube  divergence.  There  is  also  work  being  done  on  more  comprehensive 
models  allowing  for  finite  flux  values  out  to  the  equatorial  plane,  thus  des¬ 
cribing  inter-hemispheric  transport  of  ionization  (Young  et  al.,  1980a,  b). 

At  low  latitudes,  diffusive  equilibrium  is  assumed  to  be  an  adequate  des¬ 
cription  of  the  plasma  density  distribution  for  most  local  times,  and  at  mid 
latitudes  up  to  an  L-value  of  about  L*3,  diffusive  equilibrium  exists  for  a 
smaller  fraction  of  local  times;  it  can,  however,  be  regarded  as  a  limiting 
case  for  plasma  distribution  along  flux  tubes  resulting  in  a  maximum  value  for 
the  flux  tube  plasma  content.  At  higher  latitudes  it  is  not  apparent  that 
diffusive  equilibrium  ever  exists  since  the  flux  tubes  now  approach  the  plas- 
mapause  where  the  large  volume  of  the  flux  tube  combined  with  dynamic 
processes  removing  plasma  at  the  equatorial  plane  rarely  allow  time  for  diffu¬ 
sive  equilibrium  to  be  established.  Experimental  observations  related  to  the 
average  position  of  these  regions  of  the  plasmasphere  have  been  discussed  by 
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Raitt  and  Dorling  (1976). 

It  is  the  purpose  of  this  paper  to  utilize  Chiu's  formulation  of  the 
plasma  density  distribution  to  show  how  the  effects  of  flux  tube  divergence 
can  be  explicitly  allowed  for  in  terms  of  a  dipole  magnetic  field;  and  then  to 
use  this  formulation  with  experimentally  derived  boundary  conditions  to  make  a 
parametric  study  of  the  plasmaspheric  plasma  density  distribution  for  the  lim¬ 
iting  diffusive  equilibrium  case  from  L-values  of  L-1.2  to  L=3.0.  The  range 
of  boundary  conditions  show  marked  solar  cycle,  seasonal  and  diurnal  varia¬ 
tions  beyond  that  used  by  Chiu  et  al  (1979).  These  dependencies  are  difficult 
to  reconcile  with  a  diffusive  equilibrium  model,  however  by  comparing  the  Chiu 
model  with  the  results  of  the  more  classic  flux  tube  divergence  free  diffusive 
equilibrium  models  it  is  possible  to  demonstrate  that  a  diffusive  equilibrium 
formulation  is  useful  for  studying  the  variation  of  the  limiting  values  of 
total  plasmaspheric  content . 

Theoretical  Model 

The  computations  presented  in  this  study  will  be  based  on  a  hydrodynamic 
model  of  the  plasma  distribution  along  a  flux  tube  linking  the  ionosphere  to 
the  equatorial  plane.  Earlier  hydrodynamic  models  referred  to  above  did  not 
account  for  the  divergence  of  a  magnetic  flux  tube  as  it  was  followed  away 
from  the  lower  ionosphere  to  the  equatorial  plane.  In  many  cases  this  omis¬ 
sion  was  not  very  significant  since  the  plasma  distribution  was  being  followed 
only  over  a  relatively  small  altitude  range  compared  to  the  range  being  stu¬ 
died  presently.  However,  this  effect  was  allowed  for  by  Chiu  et  al,  (1979)  in 
their  hydrodynamic  formulation  by  means  of  a  force  parallel  to  B  acting  on  the 
plasma  in  the  flux  tube.  In  our  study  we  have  started  from  the  standard 
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hydrodynamic  continuity  and  momentum  equations  for  multicomponent  plasma  and 
allowed  for  the  flux  tube  divergence  by  including  its  effect  in  the  formula¬ 
tion  of  the  pressure  gradient  term  in  the  momentum  equation. 

If  we  assume  that  plasma  transport  occurs  only  parallel  to  B,  and  we 
identify  a  spatial  coordinate  S  with  distance  along  a  diverging  flux  tube,  we 
can  write  the  steady  state  equations  of  continuity  and  momentum  transport  for 
charged  particle  species  j  as: 

Continuity: 

AU)  I?  [nj  uj  A(s)1  "  Qj(s)  '  Lj(s) 


where  A(s)  «  area  of  flux  tube  at  distance  S  from  lower  boundary 


"j  «  number  density  of  species  j 

uj  *»  flow  speed  of  species  j  parallel  to  B 

Qj  *  production  rate  of  species  j 


Lj  *  loss  rate  of  species  j 
Momentum : 

S“‘  ♦  .,!■■■  [Jtf  +  g(r)  1  »s 
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where  the  additional  parameters  are: 
i® j  *  mass  of  species  j 

P j  *  kinetic  pressure  of  species  j 

q^  «  electric  charge  of  species  j 


fl£r  sin#  ec*s 


(2) 


E  *  polarization  electric  field 

g(r)  ■  gravitational  acceleration 

■  angular  speed  of  earth's  rotation 

e  *  colatitude  of  point  on  flux  tube 

r  ■  distance  from  center  of  earth 

A 

ec  *  unit  vector  of  the  centrifugal  force 

vjk  *  collision  frequency  between  species  j  &  k 

If  we  now  consider  only  diffusive  equilibrium  that  is  Uj  *  *  o  and  use  the 

fact  that  B(s)/B(so)  is  proportional  to  A(s0)/A(s)  we  obtain  the  equation  for 
the  distribution  of  plasma  density 

s 

n(s)  -  nQ  [Bs  T0/B0  Ts]'exp  ds' 

so 

tqeE*s'  +  lmQ?r  sin20  e  «s'  -  m(GM  /r3  -  sin20)r.s'j/kT(s') 

7  t  0  t  t  (3) 

directly  without  invoking  Chiu’s  parameter  F  .  Where  the  additional  parame¬ 
ters  are: 

«  geomagnetic  field  at  point  s  from  lower  boundary 
B0  «  geomagnetic  field  at  lower  boundary 
G  ■  universal  gravitational  constant 


M..  ■  mass  of  earth 


e  ■  unit  vector  of  the  transverse  component  of  the  centrifugal  force 

k  -  Boltzmann's  constant 

T  =  species  kinetic  temperature 

In  order  to  recover  the  more  classical  diffusive  equilibrium  formulation 
which  would  not  contain  flux  tube  divergence  the  magnetic  field  strength  ratio 
®S/B0  in  equation  3  would  be  set  to  unity. 

Assuming  a  dipole  magnetic  field,  the  electron  and  ion  densities  along 
the  flux  tube  can  be  given  in  terms  of  the  colatitude  (9)  of  the  point  on  the 
flux  tube,  in  the  form 

ne(s)  -  neo  [Bs  T0/Bq  Tg]  exp  <meI)/exp(4») 

(4) 


ni(s)  -  nlo  [Bs  T0/B0  Ts]  exp  (♦)  exp^I) 


where  the  integrals  I  and  *  are  defined  as 

e 

!  „  de"(-2GMccose7r0  sin39' 

9 

o 

+  3q2  r*  cose'  8inV)/kT(e') 


♦ 


A 

[qe  E  •  s'/kT(s')]  ds' 


(6) 

(7) 


Using  the  neutrality  assumption  at  each  point  along  a  flux  tube 
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(8) 


Che  integral  4  can  be  expressed  in  terms  of  I  as 
exp(4)  -  [neQ  exp(meI)/S  nlo  expfajl)]1^2 


(9) 


Combining  equations  (4),  (5)  and  (9)  it  can  be  seen  that  the  electron  and 
ion  densities  depend  only  upon  the  integral  I  the  evaluation  of  which  requires 
a  plasma  temperature  profile  along  the  flux  tube.  We  have  adopted  the  func¬ 
tional  form  of  the  temperature  used  by  Chiu  et  al,  (1979)  as 

T(s)  "  Tq  +  TL  ((L-1)/L0)a  (s/l)B 


where  l  is  the  length  of  a  flux  tube  from  lower  boundary  to  the  equator,  the 
constants  a  and  8  define  the  shape  of  the  temperature  profile  and  are  dis¬ 
cussed  later  when  actual  boundary  conditions  are  defined. 


The  integral  I  can  be  evaluated  by  Simpson's  rule  for  each  0  value  along 
a  flux  tube,  the  ion  densities  and  integral  4  can  be  developed  from  equations 
(5)  and  (9)  respectively.  In  order  to  get  the  polarization  electric  potential 
4  from  4  the  differential  relation  between  electric  field  E  and  electric 
potential  4  is  used  and  4  is  expressed  in  terms  of  4  in  the  form 
8 


4 


ds'[KT(s')d4  ] 
7b' 


(11) 


which  again  can  be  integrated  by  Simpson  rule. 


Boundary  Conditions 


The  diffusive  equilibrium  model,  described  in  the  previous  section, 
requires  density  boundary  conditions  for  all  three  ion  specie^  at  the  foot  of 
the  field  line.  In  addition  the  full  temperature  profile  from  the  foot  of  the 
field  line  to  the  equator  is  required.  For  this  study,  these  boundary  condi¬ 
tions  were  determined  empirically,  based  on  experimental  data.  The  observa¬ 
tions  of  both  ion  density  and  plasma  temperature  in  the  region  L  <  3  have 
shown  wide  dynamic  ranges  which  have  been  attributed  to  solar  cycle,  seasonal, 
diurnal  and  geomagnetic  storm  effects.  Since  it  is  the  purpose  of  this  study 
to  test  whether  diffusive  equilibrium  can  resonably  be  used  to  describe  the 
inner  plasmasphere  our  selection  of  boundary  conditions  should  reflect  the 
changes  due  to  long  period  mechanisms,  for  example,  boundary  condition  changes 
on  storm  time  scales  would  not  be  consistent  with  a  diffusive  equilibrium 
model  and  any  agreement  with  plasmaspheric  densities  would  be  purely  fortui¬ 
tous.  The  selection  of  boundary  conditions  will  therefore  concentrate  on 
reproducing  observed  solar  cycle,  seasonal  and  diurnal  variations.  It  could 
be  argued  that  diurnal  time  scales  are  also  too  short  for  diffusive  equili¬ 
brium  to  be  maintained,  however,  for  the  range  of  L-values  studied,  we  assumed 
that  both  day  and  night  ionospheres  were  in  diffusive  equilibrium.  The  prob¬ 
lems  associated  with  seasonal  effects  resulting  in  the  boundary  conditions  at 
one  foot  of  a  field  line  differing  from  its  other  foot  will  be  addressed  In  a 
later  section;  in  our  initial  study  we  assumed  that  the  plasma  distribution  on 
a  field  line  can  be  computed  for  the  segment  from  the  ionosphere  to  the  equa¬ 
torial  plane. 
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1 .  Ion  Density  at  675  km. 

The  choice  of  an  altitude  at  which  to  start  the  diffusive  equilibrium 
model  is  somewhat  arbitrary.  Since  we  intend  to  model  H+,  lit  and  0+  we  need 
to  select  an  altitude  at  which  measurements  of  all  three  ions  are  available 
over  a  range  of  solar  cycle  conditions.  An  altitude  of  675  km  was  selected, 
and  where  possible  data  was  used  only  if  it  was  within  about  25  km  of  this 
altitude.  The  relatively  small  scale  heights  for  the  heavy  ion  species  and 
the  plasma  temperature  variation  for  the  light  ions  in  this  altitude  range 
makes  extrapolation  of  densities  to  other  altitudes  prone  to  error.  Data  from 
a  variety  of  satellites  analyzed  by  a  number  of  authors  were  compared  and 
reduced  to  a  single  set  of  values  for  each  solar  cycle  -  season  -  diurnal  com¬ 
bination.  Although  this  range  of  conditions  at  675  km  is  not  completely 
covered  by  data  in  the  literature,  it  is  more  complete  than  at  other  alti¬ 
tudes.  To  help  fill  the  gaps  in  the  availability  of  boundary  conditions  and 
resolve  discrepancies  our  selection  has  been  biased  by  the  following  set  of 
currently  accepted  physical  guidelines: 

1)  Solar  minimum  densities  at  675  km  are  lower  than  at  solar  maximum. 

2)  A  seasonal  anomaly  in  F-region  plasma  density  exists. 

3)  Nighttime  densities  are  generally  lower  than  daytime  densities. 

In  searching  through  the  literature  examples  contradicting  these  guidelines 
were  found,  however,  the  bulk  of  the  data  are  consistent  with  these  trends. 
Table  1  lists  the  authors  and  satellites  used  in  obtaining  the  density  bounda¬ 
ry  conditions.  From  these  measurements  the  densities  shown  in  Table  2a  and  2b 
were  inferred.  Many  of  the  entries  in  Table  2  are  described  relative  to  one 
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particular  entry.  This  is  particlarly  true  for  composition  and  is  used  to 
indicate  trends.  At  other  times  when  no  particularly  significant  trend  was 
present,  entries  are  referred  to  as  being  the  same.  From  the  measurements  it 
was  evident  that  the  ionosphere  at  675  km  could  be  broken  into  two  regions;  a) 
from  L  ■  1.15  to  1.5  and  b)  from  1.5  to  3.0.  Beyond  L  =  3.0  effects  such  as 
the  light  ion  trough  indicate  that  the  description  of  the  plasma  distribution 
of  the  light  ions  by  diffusive  equilibrium  is  inappropriate.  In  the  inner 
region,  see  Table  2a,  the  boundary  conditions  are  latitude  independent  (L 
value)  whereas  for  the  outer  region,  see  Table  2b,  latitudinal  variations  are 
important . 

2.  Plasma  Temperature 

The  temperature  "boundary"  conditions  requires  a  full  temperature  alti¬ 
tude  profile  for  each  set  of  conditions.  Such  observations  are  not  available, 
hence  a  combination  of  ionospheric  and  plasmaspheric  observations  were  used  to 
represent  extremes  of  temperature  in  these  different  locations  and  then  a  set 
of  coefficients  for  equation  10  were  determined  such  that  a  reasonable  alti¬ 
tude  profile  of  temperature  is  obtained  using  equation  10.  The  validity  of 
equation  10  is  tested  and  discussed  in  a  later  section.  In  order  to  determine 
the  plasma  temperature  both  the  ion  and  electron  temperatures  are  required. 
Table  3  lists  the  authors  and  satellites  from  which  the  temperatures  spanning 
the  solar  cycle,  season  and  day-night  conditions  were  obtained.  These  plasma 
temperatures  and  their  latitudinal  variations  are  given  in  Tables  4a  and  4b. 
Figure  l  represents  several  cases  of  temperature-altitude  variation  for  L  *  2. 
These  profiles  are  representative  of  the  variations  contained  in  Tables  4a  and 
4b.  In  each  case  the  plasma  temperature  at  the  foot  of  the  field  line  is 
lower  than  that  at  the  equator  with  temperature  increasing  relatively  rapidly 
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at  low  altitudes.  Such  an  altitude  shape  is  consistent  with  the  limited  low 
altitude  temperature  gradient  measurements  available  and  the  calculations  from 
more  rigorous  full  flow  models  (Richards,  private  communication  1982). 

Diffusive  Equilibrium  Plasmaspheric  Densities 

The  model  described  above  was  used  to  compute  the  plasma  number  density 
distribution  along  magnetic  flux  tubes  from  675  km  to  the  equator  for  the 
field  lines  described  by  L  *  1.15,  1.5,  2.0,  2.5  and  3.0.  For  each  field  line 
a  total  of  eight  runs  consisting  of  permutations  of  solar  cycle  (max  or  min), 
seasonal  (winter  or  summer)  and  diurnal  (day  or  night)  boundary  conditions 
were  made.  These  boundary  conditions  were  determined  from  Tables  2  to  4. 
Figure  2  contrasts  two  such  runs  for  L  *  2.  Both  are  for  daytime  solar  maxi¬ 
mum  conditions  with  the  upper  being  for  summer  and  the  lower  for  winter.  Each 
panel  shows  the  H+,  He+  and  0+  density  variation  along  the  magnetic  field  line 
as  a  function  of  altitude.  Particularly  noticeable  is  the  enhanced  summer  0+ 
densities  and  the  reduced  summer  He+  densities.  These  differences  being  asso¬ 
ciated  with  both  the  differences  in  boundary  conditions,  and  with  the  enhanced 
winter  He+  being  the  winter  He+  bulge.  For  IT*"  the  boundary  conditions  are 
very  similar  in  both  winter  and  summer  and  yet  a  significantly  enhanced  summer 
profile  is  present.  This  is  a  result  of  the  interaction  between  the  three 
positive  ions  and  the  electron  fluids  through  the  ambipolar  electric  field. 
The  higher  temperature  in  summer  can  also  result  in  a  small  enhancement  of  the 
summer  profile.  Figure  2  highlights  a  major  difficulty  in  the  diffusive  equi¬ 
librium  formulation,  namely  that  the  equatorial  densities  for  the  same  field 
line  (half  summer,  half  winter)  has  different  values  depending  upon  from  which 
hemisphere  the  solution  is  derived.  This  point  will  be  addressed  in  detail  in 


a  later  section. 
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In  addition  to  the  runs  made  for  specific  solar  cycle,  seasonal  and  diur¬ 
nal  conditions,  a  study  was  made  of  the  sensitivity  of  the  model  to  boundary 
conditions  for  the  L=3  solar  minimum-summer-daytime  situation.  The  effect 
upon  equatorial  density  and  flux  tube  content  of  reasonable  variations  in 
either  the  ionospheric  density  boundary  conditions  or  temperature  profile  were 
tested.  The  equatorial  density  was  found  to  vary  almost  proportionally  to  the 
square  root  of  the  product  of  the  0+  and  H+  lower  boundary  density.  This 
strong  coupling  between  the  light  and  heavy  ions  give  a  wide  dynamic  range  of 
equatorial  densities.  In  addition  to  this  degree  of  freedom,  the  choice  of 
temperature  profile,  either  absolute  or  relative  altitude  distributions, 
yields  a  factor  of  two  increase  or  decrease  in  the  H+  flux  tube  content.  This 
factor  can  also  be  achieved  by  keeping  the  ionospheric  and  equatorial  tempera¬ 
tures  constant  and  changing  the  profile  shape.  Based  upon  these  studies  it 
appears  that  the  diffusive  equilibrium  solution  will  give  a  wide  range  of  var¬ 
iation  due  to  the  difficulty  of  supplying  sufficiently  specific  boundary  con¬ 
ditions,  and  hence  the  results  to  be  presented  should  be  viewed  as  mean  trends 
bracketing  extreme  situations  rather  than  specific  values  for  the  given  set  of 
geophysical  parameters. 

Figure  3  shows  contour  plots  of  the  H+  density  for  L  =  1.15  to  L  ■  3. 
Each  panel  shows  a  noon-midnight  meridian  cross-section  of  the  earth  and  plas- 
masphere  for  the  northern  hemisphere,  the  dayside  Is  on  the  left  side  of  each 
panel  and  the  seasonal  condition  is  indicated  by  the  shading  on  the  earth. 
Dark  shading  indicated  the  region  is  in  darkness.  The  upper  two  panels  cor¬ 
respond  to  solar  maximum  conditions  with  summer  and  winter  being  the  upper  and 
lower  of  these  two  panels  respectively.  Solar  minimum  conditions  arc  shown  in 
the  lower  two  panels  with  summer  and  winter  again  being  shown  by  the  upper  and 
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lower  of  Chese  two  panels  respectively.  For  each  plasmaspheric  condition  the 
logarithm  of  the  H+  density  along  the  L  =  1.15,  1.5,  2.0,  2.5,  3.0  field  lines 
were  combined  and  contoured  at  0.333  log1Q  (density)  intervals. 

For  the  eight  cases  shown  in  figure  3  the  largest  systematic  variation  is 
between  corresponding  solar  maximum  and  solar  minimum  conditions.  During 
solar  maximum  conditions  (upper  two  panels)  the  densities  are  between  two  and 
ten  times  larger  than  at  solar  minimum,  with  the  largest  differences  occurring 
at  the  greatest  radial  distances.  Smaller  density  variations  are  observed 
between  winter-summer  or  day-night,  in  both  instances  the  variations  are 
between  factors  of  one  and  three.  During  solar  maximum  winter  (top  panel)  the 
plasmaspheric  density  distributions  are  most  uniform,  with  under  two  orders  of 
magnitude  density  change;  in  contrast  solar  minimum  shows  over  three  orders  of 
magnitude  change.  In  general  the  H+  density  tends  to  decrease  with  radial 
distance  and  with  latitude,  although  beyond  a  radial  distance  of  about  2.5  Re 
the  density  increases  slightly  with  latitude.  These  H+  densities  are  not  sim¬ 
ply  related  to  the  H+  boundary  conditions.  As  shown  by  figure  2,  the  other 
ions,  especially  0+,  have  a  significant  effect  upon  the  equilibrium  solution. 
Indeed  the  IT*-  density  shows  proportionality  to  the  square  root  of  the  0+  boun¬ 
dary  density. 

Figure  4  shows  solar  maximum  He+  plasmaspheric  density  distributions  for 
summer  (top  panel)  and  winter  (lower  panel).  The  densities  are  contoured  at 
0.333  intervals  on  a  logarithmic  scale,  and  densities  are  labelled  in  logar¬ 
ithmic  values  of  the  density  in  cm-^.  These  two  panels  correspond  to  the  H+ 
top  panels  of  figure  3.  Comparing  the  H+  and  He+  for  summer  conditions,  top 
panel  of  figure  4  shows  that  He+  has  considerably  lower  densities,  and  a 
steeper  radial  rate  of  density  decrease  although  the  latitudinal  variations 
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follow  similar  trends.  In  marked  contrast  the  winter  He+  density,  lower  panel 
of  figure  4,  shows  increasing  value  with  increasing  latitude,  and  winter  den¬ 
sities  an  order  of  magnitude  higher  than  the  corresponding  summer  values. 
This  winter-summer  He+  variation  is  sharply  in  contrast  to  that  of  H+  and  is 
due  entirely  to  the  high  latitude  winter  helium  bulge  in  the  neutral  atmos¬ 
phere  resulting  in  the  enhanced  winter  He+  ionospheric  densities  at  high  lati¬ 
tude. 


From  both  figures  3  and  4  it  is  evident  that  not  only  is  there  a  solar 
cycle  dependence  for  the  plasmaspheric  density  distribution  but  there  are  also 
seasonal  and  day-night  dependencies.  In  order  to  study  these  latter  two 
dependencies  more  quantitatively  figure  5  shows  the  total  ion  flux  tube  con¬ 
tent  as  a  function  of  L  value  for  the  different  conditions.  The  flux  tube  is 
assumed  to  have  an  area  of  1  cm^  at  675  km  and  extends  to  the  equator,  all 
three  ions  H+,  He+  and  0+  are  included  in  the  total  density.  For  solar  maxi¬ 
mum,  left  panel  of  figure  5,  the  flux  tube  content  is  almost  independent  of  L 
value,  however  day  contents  are  larger  than  night  and  summer  contents  larger 
than  winter.  Such  a  day-night  trend  also  exists  at  solar  minimum  but  the 
solar  minimum  conditions  show  a  decrease  of  content  with  increasing  L  value. 

The  significantly  reduced  flux  tube  content  for  nighttime  conditions  com¬ 
pared  to  daytime  conditions  can,  in  part,  be  attributed  to  the  choice  of  the 
boundary  condition  altitude  at  675  km.  At  this  altitude  production  and  loss 
processes  are  still  important  for  the  light  ions  and  consequently  a  day/night 
difference  would  be  expected  which  would  not  necessarily  reflect  the  change  in 
flux  tube  content.  Indeed  the  night  time  values  must  be  treated  as  lower  lim¬ 
its  at  best  and  probably  reflect  a  boundary  condition  nighttime  density  main¬ 
tained  by  plasma  flow  from  the  plasmasphere  to  the  ionosphere.  Of  course  such 
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a  non-equilibrium  flow  situation  severely  constrains  the  validity  of  a  diffu¬ 
sive  equilibrium  formulation  should  the  diurnal  modulation  be  indeed  a  factor 
of  two  to  four  as  Indicated  by  the  day/night  differences  shown  in  figure  5. 

A  second  area  in  which  the  diffusive  equilibrium  model  is  severely  tested 
is  that  of  winter/summer  differences.  Clearly  our  calculations  which  yield 
two  sets  of  unequal  densities  at  the  equator  are  not  self-consistent.  In  fig¬ 
ure  5  it  is  evident  that  the  most  extreme  summer-winter  differences  occur  dur¬ 
ing  the  day,  the  time  when  the  largest  summer-winter  ionospheric  differences 
would  be  expected  due  to  differences  in  production.  Nighttime  differences  are 
smaller  and  should  reflect  differences  in  the  maintenance  processes  for  the 
ionosphere  which  if  controlled  by  downward  fluxes  from  the  plasmasphere  could 
be  quite  similar.  These  daytime  summer-winter  differences  would  for  solstice 
conditions  require  more  comprehensive  interhemispheric  models  to  resolve  the 
equatorial  differences.  From  the  above  discussion  we  expect  that  the  diurnal 
ebb  and  flow  of  the  flux  tube  contents  associated  with  the  diffusive  equili¬ 
brium  model  and  empirical  boundary  conditions  would  on  the  dayside  take  values 
sightly  lower  than  those  shown  in  figure  5  while  at  night  they  would  be  signi¬ 
ficantly  higher  than  that  shown  although  possibly  still  appreciably  lower  than 
during  the  day. 

Comparison  With  Flasmaspheric  Observations 

Plasmaspheric  parameters  are  observed  by  two  quite  different  techniques, 
the  first  being  the  ground-based  observation  of  whistler  dispersion  signatures 
and  the  second  being  an  in-sltu  satellite  measurement  of  a  local  parameter 
such  as  density,  temperature  or  composition.  The  whistler  technique  uses  the 
time  dispersion  of  different  frequency  waves  propagating  in  a  whistler  duct 
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(flux  tube)  from  one  Ionosphere  to  the  conjugate  ionosphere.  This  dispersion 
signature  when  combined  with  a  model  of  the  flux  tube  density  distribution 
yields  parameters  such  as  flux  tube  content  or  equatorial  density,  see  Park 
(1972)  for  a  detailed  discussion  on  this  technique.  Park's  analysis  shows  how 
a  variety  of  different  diffusive  equilibrium  models  would  produce  an  equiva¬ 
lent  variety  of  flux  tube  contents  and  equatorial  densities.  His  diffusive 
equilibrium  models  range  from  the  classical  variety  with  no  allowance  for  flux 
tube  divergence  to  "collisionless"  which  is  almost  identical  to  the  Chiu  et  al 
(1979)  and  our  own  model  which  contain  flux  tube  divergence.  The  major 
difference  between  these  diffusive  equilibrium  models  is  that  for  the  classi¬ 
cal  model  the  equatorial  density  is  aproxlmately  equal  to  the  flux  tube  con¬ 
tent/total  volume  whereas  in  the  case  of  a  "collisionless”,  divergent  flux 
tube  models  the  equatorial  density  is  less  than  20%  of  the  average  flux  tube 
content.  In-situ  techniques  of  observing  local  parameters  are  unable  to  yield 
flux  tube  content  or  average  flux  tube  density  parameters  but  instead  yield 
absolute  local  values. 

1 •  Flux  tube  content 

In  figure  5  the  range  of  model  values  for  flux  tube  content  are  shown. 
For  solar  maximum  (left  panel)  the  content  ranges  between  2  x  lO1^  to  2  x  1013 
for  all  flux  tubes  between  L  =  2  and  L  =  3.  At  solar  minimum  this  range  is 
considerably  lower,  from  6  x  10^  to  2  x  10^,  an  order  of  magnitude  lower. 
In  both  solar  cycle  cases  shown  in  figure  5  there  is  a  summer-winter  and 
day-night  variation,  with  the  day-night  variation  being  the  larger  of  the  two 
in  general.  These  model  values  can  be  compared  to  the  inferred  values  of  Park 
(1974)  and  Park  et  al  (1979)  where  whistler  measurements  have  been  analyzed 
using  the  classical  diffusive  equilibrium  model.  In  the  region  between  L  *  2 
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and  L  ■  3  Park  (1974)  obtained  flux  tube  contents  ranging  from  4  x  10^  to  1.5 
x  10  for  summer  solar  minimum  conditions  while  in  a  more  extensive  study 
Park  et  al  (1978)  obtained  values  ranging  from  4  x  10*^  to  3  x  lO*'*.  These 
values  which  cover  a  wide  range  of  solar  cycle,  seasonal  and  diurnal  condtions 
span  the  range  of  values  for  solar  maximum  shown  in  figure  5,  however  do  not 
reflect  the  lower  solar  minimum  values  shown  in  figure  5.  In  fact  the 
observed  values  should  be  increased  by  15%  to  allow  for  our  larger  flux  tube 
volume  since  our  model  is  based  at  675  km  and  not  1000  km  as  used  by  Park 
(1974)  and  Park  ®t  al  (1978).  The  whistler  observations  show  a  marked  season¬ 
al  variation  and  only  a  weak  solar  cycle  variation  based  on  the  Park  et  al 
(1978)  analysis.  These  observations  were  taken  over  the  same  time  period  the 
boundary  conditions  for  this  study  were  acquired.  We  have  been  unable  to  find 
measurements  to  agree  with  our  predicted  low  values  of  total  flux  tube  con¬ 
tents  at  solar  minimum.  Whistler  measurements  span  solar  minimum  times,  but 
consistently  show  higher  flux  tube  contents  than  we  predict.  This  suggests 
that  either  the  diffusive  equilibrium  model  is  not  applicable  or  that  there  or 
that  there  are  difficulties  in  interpreting  whistler  measurements  at  low 
L-'values  and  in  a  low  density  ionosphere. 


Equatorial  Density 


Figure  6  shows  a  comparison  of  model  (shaded  areas)  and  observation  of 
the  equatorial  density  as  a  function  of  the  L  value  (radial  distance).  For 
the  model  values  the  summer-winter  extremes  at  both  solar  maximum  and  solar 
minimum  are  used  to  represent  the  bounds  for  the  solar  cycle  shading  in  figure 
6.  The  equatorial  density  curve  shown  as  open  circles  joined  by  dashed  lines 
are  from  Park  (1974)  under  solar  minimum  conditions  using  whistler  observa- 
tions.  Park  et  al  (1978)  show  equatorial  densities  ranging  from  5  x  10  to 
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103  cm  3  at  L  ■  3.  Hence  the  whistler  equatorial  densities  span  a  range  which 
lies  above  that  inferred  from  our  model.  This  discrepancy  is  perhaps  not  too 
surprising  since  the  equatorial  density  is  highly  model  dependent  in  the 
whistler  analysis  technique  see  Park  (1972).  For  the  whistler  data  referred 
to  above  the  classical  diffusive  equilibrium  model  has  been  used  and  this 
would  then  infer  a  maximum  value  for  equatorial  density,  whereas  if  a  "colli¬ 
sionless”  model  were  used  values  would  be  decreased  by  up  to  a  factor  of  5. 
The  equatorial  density  data  shown  as  crosses  on  a  solid  line  in  figure  6 
represents  upper  limit  of  full  flux  tube  inferred  from  in-situ  observations  of 
Chappell  et  al  (1970)  and  Chappell  (1974).  This  curve  is  taken  from  the  0G0-5 
satellite  data  in  the  vicinity  of  the  1800  local  time  bulge  region  where  the 
densities  are  found  to  be  somewhat  higher  than  at  other  local  times.  This 
curve  can  be  seen  to  be  about  a  factor  of  3  higher  than  the  largest  model 
values  shown  in  figure  6.  For  L  values  greater  than  L  -  3  this  curve  is  lower 
than  the  whistler  inferred  value.  The  L-value  dependence  of  the  model  closely 
parallels  the  1/L4  decline  observed  by  Chappell  et  al  (1970),  however  the 
whistler  data  shows  a  considerably  more  gentle  fall  off.  The  final  group  of 
data  points  shown  in  figure  6  come  from  the  GE0S-1  satellite  (Johnson  1981) 
and  represent  upper  limits  on  the  cold  H+  (<leV)  equatorial  density  over  the 
1977  to  1978  period.  These  densities  were  obtained  from  two  independent  tech¬ 
niques,  direct  particle  distribution  measurements  and  measurement  of  the  local 
upper  hybrid  frequency  which  yield  plasma  frequency.  Although  these  observa¬ 
tions  were  taken  prior  to  sunspot  maximum  the  solar  EUV  fluxes  were  high, 
averaging  well  over  100  on  the  F10.7  scale.  Comparing  all  these  values  with 
the  model  clearly  indicates  an  absence  of  the  wide  dynamic  range  inferred  from 
the  model  based  on  low  altitude  boundary  conditions.  Indeed  the  range  of  var¬ 
iation  indicated  by  the  solar  maximum  shading  would,  if  increased  by  a  factor 


Page  21 


of  2  to  3  in  density,  cover  the  observations. 

3.  Diurnal  and  Seasonal  Variations 

Both  figures  5  and  6  show  that  the  model  results  indicate  a  marked  diur¬ 
nal  and  seasonal  variations.  The  seasonal  variation  is  questionable  since  the 
model  only  solves  for  half  a  flux  tube,  leaving  an  ambiguity  at  the  equator 
where  the  summer  and  winter  flux  tubes  meet.  This  clearly  indicates  that 
under  extreme  seasonal  conditions  a  flux  tube  could  not  be  in  diffusive  equi¬ 
librium  as  modelled  by  our  simple  equations.  To  date  only  whistler  observa¬ 
tions  have  been  able  to  accumulate  enough  data  to  infer  either  diurnal  or  sea¬ 
sonal  variations.  Park  et  al  (1978)  found  a  weak  diurnal  modulation  at  L  =  3 
where  noon  equatorial  densities  exceeded  midnight  densities  by  about  50%. 
This  difference  is  clearly  less  than  that  inferred  from  our  calculations.  For 
L-vaues  below  3  the  whistler  observations  also  show  a  seasonal  modulation  with 
December  values  of  equatorial  plasma  density  exceeding  those  in  June  by  50% 
(Park  et  al,  1978).  This  effect  is  not  related  to  the  model  season  variation 
which  requires  identical  density  ambiguities  at  the  equator  for  both  sol¬ 
stices.  The  fact  that  the  whistler  observations  see  a  seasonal  variation  is 
itself  indicative  that  diffusive  equilibrium  is  not  particularly  applicable. 


Summary 

We  have  carried  out  a  parametric  study  of  boundary  conditions  applicable 


to  diffusive  equilibrium  plasmaspheric  models.  These  results  were  then  com¬ 
pared  with  observed  or  inferred  values  of  the  equatorial  density  and  flux  tube 
content  to  yield  the  following  conclusions. 
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1)  The  major,  almost  an  order  of  magnitude,  difference  between  solar  maximum 
and  solar  minimum  conditions  obtained  by  the  model  based  on  empirical 
boundary  conditions  is  not  manifested  in  the  observations  of  plasmaspheric 
parameters.  A  resolution  of  this  problem  would  require  new  in-situ  meas¬ 
urements  at  solar  minimum  conditions  to  verify  that  indeed  negligible 
solar  cycle  variation  is  present  in  the  plasmasphere ;  or  more  specific 
dynamic  equilibrium  model  studies  to  understand  how  the  highly  solar  cycle 
dependent  topside  ionosphere  is  decoupled  from  the  plasmasphere. 

2)  Diffusive  equilibrium  models  which  solve  for  only  one  half  of  a  field 

line,  that  is  ionosphere  to  equator,  introduce  an  ambiguous  seasonal 
dependence.  Since  a  strong  seasonal  dependence  is  observed  in  plasmas¬ 

pheric  observations  it  is  highly  indicative  of  a  lack  of  diffusive  equili¬ 
brium  conditions  and  this  must  point  to  the  existence  of  a  dynamic  plasma 
equilibrium. 

3)  Diurnal  variations  inferred  from  the  model  are  somewhat  higher  than 

observed  primarily  due  to  the  nighttime  ionospheric  boundary  condition 
being  lower  than  that  which  would  best  be  associated  with  a  diffusive 

equilibrium  description  of  the  plasmasphere. 

4)  Under  solar  maximum  conditions  the  bulk  properties  of  the  plasmasphere  can 
be  described  by  a  diffusive  equilibrium  formulation  to  an  L-value  of  «3. 
A  formulation  using  no  magnetic  flux  tube  divergence  gives  an  upper  limit 
while  one  using  full  flux  tube  divergence  gives  a  lower  limit. 


5)  Beyond  L“3  diffusive  equilibrium  is  rarely  applicable  due  to  the  long  time 
required  to  reach  such  a  state  (many  days)  in  the  presence  of  almost  daily 
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storm  activity.  Under  the  few  conditions  when  it  would  apply  a  diffusive 
equilibrium  model  which  contains  magnetic  flux  tube  divergence  would  be 
the  most  appropriate. 

Although  many  major  advances  in  plasmaspheric  theories  have  been  made 
since  the  earliest  uses  of  diffusive  equilibrium  to  modelling  the  plasmasphere 
these  advances  have  not  resulted  in  a  "global  scale"  replacement  of  the  diffu¬ 
sive  equilibrium  plasmasphere  model.  The  more  recent  models  which  include 
time  dependent  effects  or  dynamic  equilibrium  conditions  would  use  exorbitant¬ 
ly  large  amounts  of  computer  resources  to  have  them  incorporated  into  global 
plasmasphere  descriptions.  There  are  major  requirements  for  a  global  plasmas¬ 
pheric  description,  namely;  interpretation  of  whistler  observations;  general 
inputs  for  magnetospheric  problems;  readily  available  method  to  deduce  bulk 
parameters  for  plasmaspheric  experiments;  and  a  baseline  against  which  high 
resolution  in-situ  observations  can  be  compared  in  understanding  other  mechan¬ 
isms  also  operating  in  the  plasmasphere.  Hence  we  feel  that  the  use  of  diffu¬ 
sive  equilibrium  models  although  highly  sensitive  to  boundary  conditions  can 
be  used  for  global  applications  out  to  L~3  if  they  are  used  as  indications  of 
dynamic  ranges  in  plasmaspheric  parameters  and  not  absolute  unique  values. 
Indeed  the  use  of  time  dependent,  dynamic  equilibrium  and  diffusive  equili¬ 
brium  models  in  conjunction  with  In-situ  observations  at  low  and  high  alti¬ 
tudes  (such  as  by  the  Dynamics  Explorer  satellites)  is  essential  to  understand 
the  plasmasphere  beyond  our  current  global  diffusive  equilibrium  standpoint. 
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TABLE  1 

Observations  Relevant  to  the  Ion  Density  Boundary  Conditions 


Authors 

Satellite 

Comments 

Taylor  et  al  (1968) 

0G0-2 

solar  minimum,  latitudinal, 
composition. 

Taylor  et  al  (1970) 

0C0-4 

seasonal,  light  ions. 

Brinton  et  al  (1970) 

Explorer  32 

0+  -  H+  transition  altitude 

Taylor  (1971) 

OGO-4 

Seasonal . 

Taylor  (1972) 

0G0-6 

Solar  geomagnetic  control. 

Hoffman  et  al  (1974) 

ISIS  2 

Solar  maximum,  day-night. 

Brelg  and  Hoffman  (1975) 

ISIS  2 

Solar  maximum,  seasonal, 
compos it on. 

Raitt  and  Dorling  (1976) 


ESRO-4 


Solar  minimum,  day-night. 


673ka  Ion  Density  Boundary  Conditio 
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TABLE  3 

Observations  Relevant  to  the  Plasma  Temperature  Boundary  Conditions 


Authors 

Satellite 

Comments 

Brace  et  al  (1967) 

Explorer  22 

Solar  minimum, 
seasonal,  ions 

day-night , 

Serbu  and  Maier  (1970) 

0G0-5 

Solar  maximum, 

day-night,  ions 

Bezrukikh  and  Gringauz 
(1976) 

PROCNOZ 

Plasmasphere  ions. 

Rich  et  al  (1979) 

S3-3 

Solar  minimum. 

electrons. 

Mahajon  and  Pandey  (1979) 

ISIS-1 

and  Explorer  22 

Solar  cycle,  day-night, 
electrons . 

Brace  and  Theis  (1981) 

AE-C,  ISIS-1 
and  ISIS-2 

Electrons. 
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Bailey,  G 


Banks ,  P . 


Banks ,  P . 


Banks,  P. 


Banks,  P. 


Banks,  P. 


Banks ,  P . 


Banks,  P. 


Bauer,  S. 


Bezrukikh 
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Figure  Captions 


Figure  1  Plasma  temperature  at  L  *  2.5  for  two  extreme  conditions,  solar  min¬ 
imum  night  (dashed  line)  and  solar  maximum  day  (solid  line)  as  a 
function  of  altitude  from  the  ionospheric  lower  boundary  to  the 
equator . 

Figure  2  Altitude  profiles  of  H+,  He+  and  0+  ion  density  for  solar 
maximum-summer-daytime  conditions  (top  panel)  and  for  solar 

maximum-winter  daytime  conditions  (lower  panel).  Both  sets  of  pro¬ 
files  are  for  the  L-2  magnetic  field. 

Figure  3  H+  density  contours  for  plasmaspherlc  field  lines  between  L*1.15  and 
L*«3.  Each  panel  consists  of  a  daytime  (left  half)  and  nighttime 
(right  half)  plasmaspherlc  cross-section.  The  density  contours  are 
at  intervals  of  0.333  on  a  logarithmic  scale  and  are  labelled  in 
loglO  (cm  '*)'  The  top  two  panels  are  for  solar  maximum  summer 
(upper)  and  winter  (lower)  conditions  while  the  bottom  two  panels 
are  for  solar  minimum  summer  (upper)  and  winter  (lower)  conditions. 

Figure  A  He+  density  contours  for  plasmaspherlc  field  lines  between  L=1.15 
and  L**3.  The  upper  panel  being  for  solar  maximum  summer  conditions 
while  the  lower  is  for  solar  maximum  winter  conditions.  The  format 
and  contour  intervals  are  the  same  as  for  figure  3. 

I 

Figure  5  Total  ion  flux  tube  content  from  the  ionosphere  to  the  equator  for 
solar  maximum  (left  panel)  and  solar  minimum  (right  panel)  as  a 
function  of  the  dipole  L  value.  In  each  panel  the  seasonal  summer 
(S)  and  winter  (U)  variations  are  shown  for  both  daytime  and  night. 
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This  seasonal  difference  is  highlighted  by  the  crosshatching. 

Figure  6  Equatorial  density  comparison  between  model  solar  cycle  conditions 
and  observations  over  the  plasmasphere  as  a  function  of  radial  dis¬ 
tance  . 
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